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Summary 

Several energy- momentum "tensors" of gravitational field are considered and 
compared in the lowest approximation. Each of them together with energy- 
momentum tensor of point-like particles satisfies the conservation laws when 
equations of motion of particles are the same as in general relativity. It is 
shown that in Newtonian approximation the considered tensors differ one from 
the other in the way their energy density is distributed between energy density 
of interaction (nonzero only at locations of particles) and energy density of 
gravitational field. 

Starting from Lorentz invariance the Lagrangians for spin-2, mass-0 field 
are considered. They differ only by divergences. From these Lagrangians by 
Belinfante-Rosenfeld procedure the energy-momentum tensors are build. Us- 
ing each of these tensors in 3-graviton vertex we obtain the corresponding 
metric of a Newtonian center in G 2 approximation. Only one of these "field- 
theretical" tensors (namely the half sum of Thirring tensor and tensor obtained 
from Lagrangian given by Misner, Thorne and Wheeler) leads to correct value 
of the perihelion shift. This tensor does not coincide with Weinberg's one (di- 
rectly obtainable from Einstein equation) and gives metric of a spherical body 
differing (in space part of metric in the first nonlinear approximation) from 
Schwarzschild field in harmonic coordinates. As a result a relativistic particle 
in such field must move note according general relativity prescriptions. 

This approach puts the gravitational energy-momentum tensor on the same 
footing as any other energy-momentum tensor. 

1 Introduction 



General relativity is a complete, elegant, and self-consistent theory. Yet there 
is a necessity to obtain gravity by field-theoretical means starting from flat 
spacetime, see e.g. [1-4]. It is widely believed that on this way even dropping 
the general covariance requirement we naturally get general relativity. It is 
supposed that in the lowest nonlinear approximation this is demonstrated in 
detail by Thirring [2]. Yet this conclusion can not be drawn from [2], see Sec. 2. 

The energy-momentum tensor of material fields in general relativity is ob- 
tained from that without gravitational field by equivalence principal (comma 
goes to semicolon) . This means that general covariance dictates the form of ver- 
tices containing material fields. Even in this case other considerations may lead 
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to modifications. So conformal invariance leads to Chernikov-Tagirov energy- 
momentum tensor [5] . Dropping general covariance gives more freedom in choos- 
ing vertices. 

Since the gravitational collapse is considered as the greatest crisis in physics 
[6] the research into possible alternative theories acquire especial significance. 
It is quite natural to make the first step and to consider the simplest processes 
by utilizing vertices; the graviton propagator is known by analogy with electro- 
dynamics. 

In the lowest nonlinear approximation it is necessary to know only 3-graviton 
vertex. We assume the simplest possibility: the source of graviton is the energy- 
momentum tensor of two other gravitons. In higher approximations probably 
other vertices will be needed. Along this path one can find out what theories 
are possible without assuming general covariance and a priori restriction on 
vertices. An important step in this direction was made by Thirring [1-2]. We 
continue his investigation in the same approximation and restrict ourselves to 
point-like classical particles as sources of gravitation. Mainly we are interested 
in the simplest system consisting of a Newtonian center and test particle moving 
in its field. 

In general relativity classical particles move along geodesies in Riemannian 
space. This is the incarnation of equivalence principle and it is more reliable 
than specific equation determining the gravitational field [9] . As to the equation 
determining gravitational field, it is possible to think that the phcnomcnological 
field-theoretical approach will lead to more complicated algorithm for getting the 
field. An interesting possibility in this direction was pointed out by Schwinger 
[10]. 

It is reasonable to believe together with Einstein that for some reason or 
other the singular behaviour near the gravitational radius does not correspond 
to reality see §15 in .Pais's book [11] and Einstein paper [12]. . At present 
the Schwarzschild singularity is considered as fictitious by many researchers 
because the geometry is nonsingular there. See however the text after (2.2.6) 
in [13] and after eq.(9.40) in [14], where they say convincingly about physical 
singularity. By field-theoretical approach it is difficult to understand why in a 
finite system the acceleration of a freely falling particle becomes unlimited when 
it nears the horizon. Such a behaviour should be connected with the fact that 
according to [15] the gravitational energy in vacuum outside the sphere of radius 
R goes to — oo for R — > r g . In conformity with this the energy of matter and 
gravitational field inside the sphere of radius R goes to +oo, in such away that 
total energy of spherical body is equal to its mass. But if a theory predicts that 
the absolute value of field energy outside sphere of radius R might be greater 
than total energy of a body then the analogy with electrodynamics suggests 
that the concept of external field becomes inapplicable [16]. The belief in general 
relativity in similar circumstances is based upon the concept of nonlocalazability 
of gravitational energy, see, e.g. §20.4 in [6]. What is more, general relativity 
does not need as a rule the gravitational energy -momentum pseudotensor. 
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The situation changes drastically, when we begin to construct gravity theory 
starting from flat spacetime and assume that in 3-graviton vertex each graviton 
interacts with energy-momentum tensor formed by two other graviton. Then 
the nonlinear correction to the motion of a test particle depends on the chosen 
energy-momentum tensor. The latter is build from field Lagrangian, which is 
not unique as one can add to it some divergence terms. This leads to different 
energy-momentum tensors. They can give rise to gravitational energy densities, 
which may have even different signs. The question of sign of energy density is of 
interest by itself. Provided the sign turns out to be positive, one should expect 
the weakening of gravitational interaction at r ~ r g = 2GM in comparison 
with Newtonian one in order that the gravitational energy outside the sphere 
of radius r were much less than the mass M of the center. The possibility of 
decreasing of interaction at small distances is suggested also by the behaviour 
of attraction force between two bodies supported by Weyl's strut, see, §35 [17]. 

In order to understand in what way the various energy-momentum tensors 
differ one from the other we consider the following tensors: Thirring's [1,2], 
Landau-Lifshitz's [16], Papapetrou- Weinberg's [9] and tensor obtained from the 
Lagrangian given in Exercise 7.3 in [6]. The second and third tensors are rep- 
resentatives of general relativity, the rest are build from Lagrangians of free 
field of spin-2, mass-0 particles and symmetrized by the Belinfante-Rosenfcld 
procedure, see, §1 Ch.7 in [18]. 

In the considered approach it is suitable to subdivide the 3-graviton vertex 
into three vertices in accordance with three possibilities for choosing two gravi- 
tons out of three to form the energy-momentum tensor,- the source for the third 
graviton, see Sec. 2. So three diagrams contribute to nonlinear correction to the 
field. In contrast with this the energy-momentum tensor, figuring in the solution 
of Einstein equation by iteration procedure, is so defined that the correction to 
field is given simply by means of propagator, i.e. by one diagram only. 

The main result of the paper is this: starting from quadratic Lagrangians 
(differing by divergence terms) of spin-2, mass-0 particles, the energy-momentum 
tensors are constructed by Belinfante-Rosenfcld procedure. It turns out that 
only certain combination of these tensors (used in 3-graviton vertex) is fitted 
for correct description of perihelion shift. This combination docs not coincide 
with Papapetrou- Weinberg tensor. 

The investigation of possibilities of phenomenological approach to gravita- 
tion without use of general covariance seems to us very promising. Valuable 
undertaking in this direction was made in [19]. 

2 Thirring's energy-momentum tensor 

Throughout the paper we use 

9!iu = Viiv + h ^ (!) 



3 



In Sections 2, 3, 5 rj^ u = diag(l, — 1, — 1, — 1, ); in Sections 4, 6 and Appedix 
rj^ = diag(—l, 1, 1, 1, ). In this Sec. we use Thirring's notation [2]; both greek 
and latin indices run from to 3. 

The gravitational field is described by the symmetric tensor h^ v , which con- 
tains spin-2 and lower spins, see, e.g. [3]. The unnecessary spins (spin-1 and 
one of spin-0) are excluded by Hilbert gauge : 

hr, v = {br-\irh\ v = ^ h = h\, h v = ^h. (2) 

One way to obtain Thirring's tensor is to start from general relativity Lagrangian 
^P^gR. If we remove terms with second derivatives of g^ v into divergence terms 
and drop the latter, we get the function G(x) in eq.(93.1) in [16]. Retaining in 
it only quadratic in terms we get 

G(x) = i[V,A^' A - 2V,A^'" + 2^/^" - h x h> x ] (3) 
This is equivalent to Thirring's Lagrangian [2] 

L= lbP^,xr v ' X - + - M' X \- (4) 

Here 

V = -W 2 /. ^ = / 2 = 8^G, G = 6.67 • 10~ 8 cm 3 /g ■ sec 2 - (5) 

Using tp^ instead of h^ v is justified because then the analogy with electrody- 
namics becomes more close: ip^ is analogous to vector- potential and has the 
same dimensionality, M\[G has the dimensionality of electromagnetic charge. 
We note that the Lagrangian (4) exactly corresponds to Schwinger's Lagrangian 
[10], who uses the notation rj^ v = diag(—l, 1, 1, 1, ) and 2h^ h — —h^ v — —h^ u . 

The canonical energy-momentum tensor following from (4), has the form 

T~* s = ip^' 5 ^ - ^<V 7 - 2^' W - rf S L 

L= \[^,x^"' X - \v,x^ X - 2<p^xV Xv >^ (6) 

^=^=V^-^^V' ; ■0 = '0<t' T - (7) 

Using ip^ instead of ip^ v is handy as many expressions become more compact 
and the consequences of imposition of Hilbert gauge more clear. 

The energy-momentum tensor for a static point-like mass (Newtonian center) 

M 

T„„= M5(x)5„ 5„o- (8) 
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In linear approximation this source generate the field 

f M 

</V — —hp„/2f = - SftoSvo, hp V — AffinoSvo, (9) 

47T|X| 

satisfying Hilbcrt condition (2). For one Newtonian center <f> = —GM/r. For 
several centers 



In terms of 

h tlv = h tlv - -rinJi, h = h a a , (11) 

we have 

V = 2 ^<W> h = K° = -40 = -ft. (12) 
The energy density of field (9) is positive 

T°° = ^(V0) 2 . (13) 

/ 

But T 7 " 5 ought to be supplemented to symmetric one by the spin part: 

0~i 5 =T~i 5 +T~< 6 . (14) 
For Newtonian center Thirring obtains 

-^(v# 2 M«o, ^-^W) 2 ^. (15) 

So in this case 00 is negative 

= ~8^G (W)2 - (16) 

Turning now to conservation laws of total energy- momentum we remind first 

how matters stand in general relativity. There the energy-momentum tensor of 

p 

point-like particles T^ v is connected with its counterpart in special relativity 
T» v by the relation, sec (33.4), (33.5) and (106.4) in [16]: 

P fl 

j— gT ^ =r ^ = Y^m a uV—5{x- x a {t)), u^^dx^/ds, (17) 

a 

g is determinant of g^ v . In terms of the conservation laws are (see (96.1) 
in [16]) 

T"„,„ = [T^inrv + h TV )l» = \h^JT»°. (18) 
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We shall see below that T^ T h Tl/ can be interpreted as (part of) interaction 
energy-momentum tensor. 

As is known the equation of motion of particles in general relativity 

^ + rV"^ = o (19) 

is contained in conservation laws. Indeed from 

T jk ;j =T jk j + T 3 mj T mk + T k mj T jm = (20) 

v ., 

taking into account that from definition of T 3 in (17) 

t jk j = -\{-grk-g)A 3k + (-9)-^" d , it . = ±. g>m , (21) 

we get 

T 3k j+Y k mj T mj =Q. (22) 
This is equivalent to (19), because [9] 

T> m j = J2^6(x-x a (t)). (23) 

a 

Going back to field-theoretical approach, we rewrite the equation of motion of 
particles (19) in the lowest approximation 

^- = ^ = h a ^u a uP ~ W a . f3 u a uP. (24) 
as ds A 2 

Just at such movement of particles the divergence of total energy-momentum 
ought to be zero, and inversely, from zero divergence follows cq. (24). From 
(23) and (24) we find 

T< s „ = \h a ^T a& ~ h s aa T^ , (25) 

This agrees with (22) and (18) in considered approximation. With the same 
accuracy this can be rewritten as 

(T7 , + jTra^^ = I/j q /T^. (26) 

Using linearized Einstein equation for ip^ 

T ltv =T liV -^T , T = %° , (27) 
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we get 

T j5 ,-y - f^' 5 f aP = ~h a ^ s T a ^ = -h a ^T a P. (28) 

Addind (26) and (28) we see that the total energy-momentum tensor is con- 
served. Spin part of energy-momentum tensor is conserved by itself and do not 
contribute to conservation laws : 

T 7<5 , 7 = 0. (29) 

So from (26) and (28) it follows that the conserved tensor contains in itself 
the interaction tensor [2] 

j-'t* =T-< a h a s . (30) 

But it is not symmetric. In order to understand the reason we have to 

consider the properties of T jS in detail despite the fact that it does not take 
part in conservation laws written in the form of eqs.(26) and (28). According 
to known rules [2,18] we have 

f jk = -pi ik ti - F kij A - F ik \i , (31) 

81 f 

F M = (<p k a r, i fl -<p i a ri k f l ), L=L. (32) 

0(Paf3,j 

The antisymmetric part of (31) is contained only in the last term. For it we 
have 

\Ofaj.iJ A \OiPak,iJ ti 0(fi ajtl Oif a k,i 

The first two terms on the right hand side symmetrize the interaction tensor, 
the last two terms symmetrize the canonical one. 

It is not seen directly from (6) and (31) that field energy-momentum tensor 
9 lS in (14) is symmetric. This agrees with the fact that the proof of symmetry 
utilizes the Euler-Lagrange equations for field which is considered as free [18]. 
We are interested in interacting field. So using linearized Einstein equation (27) 
with source, we get 



dip a j,i J A \dtp a k,i 

= f{T a H k a - T°Va) = \<J ak b? a - T^h k a ). (34) 

In the last equation we use the connection between ip^ v and h^ v , see (5). Substi- 
tuting in (30) 7 — > j,S — > k, we see that the sum of (30) and (34) is symmetric. 
This result retains if we start from another Lagrangian differing from Thirring's 
one in (6) by divergence because the linearized equation remains the same. 
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' b 

One should take into account however that the symmetric part of T 3 can 
also contain terms of interaction type. So for the Lagrangian in (6) similarly to 
(34) we find 

_ pjik . _ F kij. = [fT _ 2w «]<fi + f[T a ^ k a + T Q V'J 

-2(^' fc '<V Q + ^ k ' a v\,i) + 2^ fei 'V«,*- (35) 

s 

As a result we get for T 1 

T jk = 2[(^ a ^ k a + ^ k ' a i^ a ) - V ia , ia <p kj + V«,i - <P Ac V« 

-^ k > a y\,i + ¥> W, Va,i] ~ 2^Va, 4 + 2/T'Va. (36) 

/ .. 

Here last but one term, added to T 3 , makes it symmetric. The last term can 
be rewritten in terms of h^ u in the form, see (9) and (11), (30), 

-T 3a (h a k - ^ Va k h) = - T jk + \r ]k h. (37) 

b 1 b 

So the symmetrization of T in (30) is reduced to its replacement by ^T° k h. 
This tensor is nonzero only where particles are present. For Newtonian centers 
the corresponding energy density 

X -r m h = -2T°> (38) 

int 

is positive (contrary to our intuition and) contrary to T in (30), sec (12) and 
(10), where h and <f) are given for Newtonian centers. 

We note that the use of linearized Einstein eq. (27) in the expression for T 1 
leads to that eq. (29) is satisfied only with considered accuracy. The presence 
of interaction energy-momentum tensor means the appearance of such vertex: 
the energy-momentum tensor of matter together with one of gravitons serves as 
a source for other graviton, see Fig.l 

Now we note that Belinfante-Rosenfeld procedure leads to the appearance 
in gravitational energy-momentum tensor terms with second derivatives. 

Thirring assumes that his tensor 9 l5 (see (14), (6), (31)) is an analog of 
energy-momentum tensor figuring in the r.h.s. of Einstein equation when it- 
eration procedure is used. In other words the nonlinear correction to field is 
given only by diagram (2a) in Fig. 2. On this Fig. the short straight line 
has only conditional meaning: it represents the source of gravitons, namely 
the energy-momentum tensor build of two gravitons (real or virtual) shown on 
Fig. 2 as joining the ends of this line. The graviton emerging from the middle of 
straight line is emitted or absorbed by this source. On diagram (2a) the energy- 
momentum tensor is build from gravitons of Newtonian center. On diagrams 
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(2b) and (2c) one of the virtual gravitons of Newtonian center interact with 
energy-momentum tensor of two other gravitons. All three diagrams on Fig. 2 
correspond to one Feynman diagram obtained by contracting the short straight 
line to a point. 

The contribution to nonlinear correction for field from diagram (2a) is easy 
to obtain. Indeed, from (14), (6) and (15) we have 

^=T 3 ' fc = ^(^,k--%(V0) 2 ), j, A = 1,2, 3. (39) 

Using now the field equation in Hilbert gauge with 9^ from (16) and (39) 

d 2 

□ftf" = -167^0"", □ = — - V 2 (40) 

ut 

we find 

h m = -14>\ W k = - < ±^x i x\ <f>=- — , i,jfc = l,2,3. (41) 
Here easily verifiable relations 

2 x^_25ik 4x'x fc X7 2]__2l (ao\ 

were used. The obtained satisfies Hilbert condition. Going over to = 
h^ u — we find the following nonlinear corrections 

/ llfc = _G 2 M 2 (^ + ^). (43) 

(2) 

Index 2 in h indicating the order of correction in powers of G, is dropped 
for brevity. 

Finally from (12) and (43) we have 

ds 2 = g 00 dt 2 -(l-2<f> + Z<t> 2 )5 lk dx l dx k - <f> 2 — , (44) 

3 oo = l + 2^-40 2 . (45) 
The transition to spherical coordinates is given by the relations 

6 ik dx i dx k = dr 2 + r 2 (d9 2 + sin 2 9dip 2 ), = ( dr f. 

Why the nonlinear correction in (45) turns out to be negative? It will appear 
later on that it is caused solely by the source (15), see eqs. (95) and (93a). The 
sources (8) and (15) have different signs, but the corresponding fields have the 
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same sign. The answer is simple. The correction in (45) is only a small (and 
negative) part. The larger and positive part goes for converting initially bare 
mass in Newtonian potential into a dressed one, see eq. (A9). Now the negative 
sign of correction is clear: the mass of Newtonian center at infinitely large 
distance appears as M, but at finite distance the test particle feels a greater 
mass and greater attraction, because (15) is negative. 

The nonlinear correction — 40 2 in g 00 in (45) is of special interest for us. The 
correct value necessary to explain the perihelion shift is +2(f> 2 . The shortest way 
to see this is to use the method described in §101 in [16]. We write in spherical 
coordinates 

ds 2 = A(r)dt 2 - B(r)dr 2 - C(r)(d6 2 + sin 2 9d<fi 2 ). 
The solution to Hamilton-Jacobi equation has the form 

S = -£t + + S r (r), S r (r) - J B(r)[-^ -^-^ dr - 

Here £ and J are constants. 

For nonrelativistic particle £ = m + £' , £' <C to, and the main terms in 
square bracket in the expression for S r are cancelled out: 

-l«-20(l-0), 



A(r) 



where we have assumed that A(r) = 1 + 20 + 2</> 2 . 

So we have to retain in A _1 (r) terms of order </> 2 , but in B(r) and r 2 C(r) _1 
only terms of order <j>: B(r) = r~ 2 C(r) = 1 — 2(f). Thus 



B(r)(A(r)- x - 1) w -(1 - 20)20(1 - 0) w -20(1 - 30) = ^ 1 ^ 



2 

g_ 

2r 2 ' 



The leading term cx r in S r is 



1 3m 2 r 2 
~r 2[J 2 



As explained in [16], the expansion 



S r{ J 2 - 3 ^) = S r (J 2 )-^ dS ' 



2 d,P 

directly leads to correct perihelion shift. 

In general the motion of a particle is described by equations (cf. §4 in Ch.8 
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For nonrclativistic particle only the third term on the l.h.s. of second equation 
in (45a) does not contain small factor of order v 2 . So it requires more accu- 
rate approximation. For relativistic particle A(r), B(r) and ought to be 
considered in the same approximation. 

Now taking into account all 3 diagrams of Fig. 2 we get instead of (44) (see 
Appendix for more details and pay attention to difference in metric signature 
there) 

ds 2 = (1 + 2<p + 4<p 2 )dt 2 - (1 - 2<f) + 9<p 2 )(dx) 2 - 130 2 ^-dxW, (456) 

As the nonlinear correction to 500 is twice as much as necessary, Thirring tensor 
alone is insufficient. 

We note here that Thirring obtained from his tensor the necessary correction. 
Yet his result is objectionable as he used illdefined gauge 

2 . G 2 M 2 

see eq. (83) in [2]. Namely the source of A fall of too slowly for large r and the 
integral defining A , see (A8), diverges for large r'. 

In the next two Sections we shall see how energy- momentum "tensors" of 
general relativity differ from Thirring's tensor. 

3 Landau-Lifshitz pseudotensor of energy- momentum 

This pseudotensor in the sense and approximation considered here is a tensor. 
In the lowest approximation with help of relation 

V~99 lk « (1 + \)W k - h ik ) « r? k - h ik . (46) 
we get from eq. (96.9) in [16] 

xik _' [ZikjlZ ,m Til I km _ T km,p T ,i Tirn.p T ,k _■_ Zim,pZk , 

L ^ ^ [11 'Hm 'v ,m '^mp '*> '^mp ~r ft ft m,p > 

h^% q - k - h^h* + if k i}-h mn ^ m - n + h' m ~h m - \h m , m h^ m )]- (47) 

Comparison with canonical tensor (6) shows that it is connected with t lk by the 
relation 

t ik =T ik +F ik , h ik = -2f<p ik , 
F ik = ^-^[h ik '% n > n - h l \ih kn . n - h kn '"h n / + h in *h k n , p ). (48) 
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From (14) we see that now in place of T lk stands F lk . But T lk was a con- 
served quantity, see. (29). So F lk should rather play the role of interaction 
energy-momentum tensor. Indeed, taking into account that in considered ap- 
proximation hik satisfies the linearized Einstein equation 

h*np,j^ hjp,n^ hjn.p^ ~i~ TJnphqr ^ 167rGr^7f i p, (^9) 

we find 

F j k j = ~ h kn, Tm = h kn,i Tn . _ \ h ^r k . (50) 

Now we check that coservation laws [16] 

^((-9)[T^ k +f k })^0 (51) 
arc fulfilled. From (48), (28) and (50) we have 

t ik ,i = -\ti^ k % q + h kn ^T m - l -h l r k . (52) 
For matter energy-momentum tensor from (17) we get 

(-g)T lk = V=gT lk *(l + ^)T lk , -g*l + h. (53) 

From here 

((-g)T lk ) ^r k ^ + h t r k , (54) 

the terms of order h 2 being dropped. Now it follows from (25), (52) and (54) that 
(51) is fulfilled. As seen from (53) here too there is a tensor, which is nonzero 
only where particles are located. Surprisingly it coincides with Thirring's inter- 
action tensor, see (37) and text below it. 

Although —g T lk +t tk differs from Thirring's T tk + T lk +0 tk } for Newtonian 

centers they coincide, see eqs. (16) and (39). 

Now we turn to Newtonian approximation. According to Problem 1 in §106 

in [16] the energy density of gravitational field in this approximation is given 
/ 

by T 00 in (13), (10). But there is also energy density of interaction fi<j), where 
fi is density of particles. Using Poisson equation V 2 = AirGfi and ignoring the 
problems connected with point-like nature of particles, we can write (utilizing 
integration by parts) 

J v4W = -f i^( v< ^) 2 ^ ( 55 ) 

The density in the integrand on the r.h.s. contains now not only the energy 
density of interaction, but also the proper energy density of particle's self-field. 
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The density on the l.h.s. is nonzero only at particle locations, the density on the 
r.h.s. is nonzero where the field is nonzero. The integration by parts deprive us 
the possibility to retain the previous physical meaning of integrand. If never- 
theless we do this, then adding to (13) the energy density in the r.h.s. of (55) 
we get the effective gravitational energy density in Newtonian approximation 

1 ;(V0) 2 . (56) 



8ttG v 

To bring this in agreement with t 00 we ought, according to a foot-note in [16], 

p 

take into account the contribution from (— g) T 00 ■ Let us do it. For t 00 we have 

ta ° = -^G m)2 > (57) 
where (j) is the potential of Newtonian centers. Now 

(-.9) f 00 = V—gT 00 « r 00 (i + « r 00 - 2<m, (58) 

see (12). The sought for agreement will be reached only after we rewrite a la 
Thirring [2] T 00 in terms of observables. From (17) 

T 00 =£m„=U(af -*„(*)). (59) 

a 

In the presence of gravitational field 

ds 2 = g oa dt 2 (l-v 2 ). (60) 
Here v 2 is physical velocity, see §88 in [16]. Hence 

dx° 1 1 hoo 



v 2 < 1. (61) 



Thus after going over to the observable velocity we obtain the term 

-^% = -<M- (62) 

detached from T 00 . Equation (12) was used to get the r.h.s.. Together with 
corresponding term in (58) this leads to 

-3/^=i>-^(V# 2 , (63) 

where arrow corresponds to going over in (55) from integrand on the l.h.s. to 
the integrand on r.h.s.. Now the sum of (57) and (63) gives the expected (56). 
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The consideration of Newtonian approximation makes the following point of 
view very enticing: The energy density of an isolated point-like particle should 
be positive; Hilbert gauge exclude the unnecessary spins and then positivity 
seems quite natural, because the presence of virtual gravitons should not make 
the energy density negative. The attraction is described by interaction energy 
density and so the latter must be negative. Neither Thirring tensor nor Landau- 
Lifshitz one satisfies this requirement. The MTW tensor does. Unfortunately I 
failed to fit this idea into existing approach to gravitation. 

Using LL tesor in 3-graviton vertex we get from diagram (2a) the same 
contribution as in the case of Thirring tensor. The contribution from all 3 
diagrams of Fig. 2 leads to 

ds 2 = (1 + 2(j) + Acj) 2 )dt 2 - (1 - 2<f> + 7cf) 2 )(dx) 2 + l^^-dx'dx 3 . (63a) 

4 Papapetrou- Weinberg energy- momentum ten- 
sor. 

Einstein equation can be recast in such a way that gravitational energy -momentum 
"tensor" can be easily identified in coordinate system that goes over to Minkowski 
system at large distances from gravitating bodies [9] . In the lowest approxima- 
tion this tensor have the form, see eq. (7.6.14) in [9]) 

(64) 

In this Section we use Weinberg notation: 
r^,, = diag(-l, 1, 1, 1), = ^ v + h^, = h™ v = -h Thtr = -h^. (65) 

Greek indices run from to 3, latin ones from 1 to 3. R^i) 2 ^ is Ricci tensor in 
the first and second approximation in powers of h^ v . The indices are raised and 
lowered with -q^. In terms of we get 



1__ _ _ _ 1__ 

^ (^^K^ ^^^cr,M ^^M^^) 2^ tf^ih fi.K "i" h n ^ 
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+^ K [\h Xv h M - ±h~h x x + X -h v a , v h'° - ±h,M], (67) 



+\h Ka ,xh KX >° - \h Ka ,xh K °< x - \hh KX , KX + \h.M- (68) 



For Newtonian center from (9-12) we obtain 

>4S^, h = h w = h T = -40 = 

(69) 



hpv = -h£ v = -40^ o <Lo, h^ u = -ti^ v = -2(j)S^, h = h w = h T = -4(f) = -h. 



Nonzero components of t flK are 

3 ,_,, 2 3GM 2 m GM 2 2 ^ . k . 
tm = _ 8^G (W) = ~^r7^' Uk = 8^ (?<W _ UX X ) - (70) 

In Hilbert gauge from equation 

V 2 V = -16ttGV (71) 

we find, cf. with (40), (42), 

1G 2 M 2 T i r k 
~ h ™ - ^F- - 30 2 , ji* = _7G 2 M 2 ^. (72) 

It is easy to check that (72) satisfies the Hilbert condition (2). In terms of 
we have 

h 00 = -20 2 , h ik = G 2 M 2 (^ - 7^p). (73) 

In the expressions (71-73) h^ v is nonlinear correction. 

On the other hand in harmonic coordinates the Schwarzschild solution has 
the form [9] 



-dr 2 = -Y^dt 2 + (1 - 0) 2 (df) 2 + t-^^ 2 ^ dx ^ dx \ (74) 

So in the considered approximation this gives 

.goo = -(l + 2^ + 20 2 ), (75) 

^ = (1-2^ + ^ + ^). (76) 

From (69) we have /iqq = —20, from (73) h$ — — 20 2 , and there is agreement 
with (75). As to the nonlinear correction for gik, in (73) it differs from the one 



15 



in (76) by a gauge. Really, subtracting from hit in (73) the nonlinear part of 
(76), we find 

G 2 M 2 (^ - ^) = 2G 2 M 2 (A, fc + A M ), A< = ^, (77) 

i.e. a gauge. 

Going back to too in (70), we note that this density is negative and does not 
coincide with any density of other tensors. At the same time the equation of mo- 
tion of particles is contained in the conservation laws of total energy- momentum 
tensor. We shall check it in considered approximation. For gravitational part 
the calculations give 

^% = -h\^ a + h a T^ - h?^T pa - h vX T» ViX . (78) 

The energy-momentum tensor for particles, figuring in conservation laws, has a 
rather complicated form by construction [9] 

r"" = rfW^r/j T a(3 ~T ^ K + W a T aK + h K a T ap 

r, T" K - hiT^ K + h" a T aK + h K a T a ^. (79) 

From here with the considered accuracy 

t^, k =T ^ K , K + W a , K T aK + h K a . K T a " + h K a T a ^ K . (80) 

So 

(r^ + t» K ) tK =T »\ K - \h pa ^ pa + h^- K T aK + h a T^. (81) 
Further from (53) we get 

T « T^, K - \h K T» K . (82) 

Taking into account (25) we see that the r.h.s. of (81) is zero. 

Now we turn to Newtonian approximation. Terms of interaction tensor are 
contained in both t^ v (three last terms in the r.h.s. of (79)) and in t^. From 
(64), (66-68), using (49), which preserve its form in the notation of this Section, 
we find the following terms of interaction tensor contained in t pK : 

-\h pK T - ^ K Qi pa r pa - -~hT) - ^hT pK . 

From here and the first equation in (70) we get in Newtonian approximation 

£oo = -^(V0) 2 -60T O o. (83) 
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Here <j) is the same as in (10). From (79) and (53) we get in this approximation 

T oo = (-g)-iT 00 + 2h%T m w T 00 (l - h) - 2h m T m = T 00 + 60T 00 . (84) 

Thus in Newtonian approximation the interaction terms in the sum of (83) and 
(84) are cancelled out. The agreement with Newtonian approximation (56) is 
achieved in the same way as for Landau-Lifshitz tensor: T 00 on the r.h.s. of (84) 
detaches term (62), which is equivalent (in accordance with (63)) to 4^j(V0) 2 . 
Together with the first term on the r.h.s. of (83) this gives (56). 

Weinberg shows in detail that his energy-momentum tensor has all required 
characteristics. But this tensor does not help us to find energy-momentum 
tensor of two gravitons as represented by straight line on diagrams of Fig. 2. By 
construction Weinberg's tensor gives the gravitational field only via diagram 
(2a). The field-theoretical description tell us that test particle is not quite 
passive. It does not simply follows the command "move along geodesic" but 
itself takes part in the creation of field in which it moves, see Fig. (2b, 2c). 
From this viewpoint one can expect that e.g. photon and graviton scatter 
differently on Newtonian center as only in the latter case all three diagrams of 
Fig. 2 contribute. 

5 Misner-Thorne- Wheeler energy- momentum ten- 
sor 

In this Section it is handy for us to use again Thirring's notation. Up to diver- 
gence terms the Lagrangian (6) may be rewritten as [6] 

C = \p^\^ X - - W'V",*- (85) 

The corresponding canonical energy-momentum tensor 

T jk =g^r-^ v ' k -rf k C, (86) 

acquires the form 

T jk = f jk -^rf k f,f ik = <p^ k tp^ J 2<{P v > k <p va ' a ,T = —T = 2C. 

(87) 

Spin part is given by (31-32) with substitution L — > £. We dwell on differences 
from Thirring's tensor. In symmetric in jk tensor 

F jik + pktj = ^ai,j _ ip oca ^ lf) k a + ^ai,k _ ^ ^ki^j q 

-2^%^' + {2^\„rf k - y ak ' j - Va + V ka ,*V ij + ^'%¥> ifc (88) 
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there is no derivatives over x % . This means that in divergence -F- 7 ,j + F kl \i 
there are no terms of interaction tensor. In antisymmetric in jk tensor 

Fikj / ak,i aa j ka ij , / aa ij aj.i\ k , ja ik /oc\\ 

J — [IP ' ~<P ,aV )ip J a — <P ,aP +{¥ ,<jf\ -<P )<P a+p J ,aP (89) 

such terms are present. Hence, using linearized Einstein equation (27) in the 
expression for F %k i ^ , we get 

_F*i. = /(f - T fe V' Q ) + if" a * - y k J). (90) 

Terms with / together with (30) give symmetric interaction tensor in accordance 

/ 

with (34). Other two terms on the r.h.s. of (90) supplement T jfe to symmetric 
one, see (87). So we get 

Qjk = jk + j, jk = ^k^j _ l^kpj _ ^>-^k _ pkv ^ij 

-¥> Qi Va,i - fV«,i + {V> ai,k i + V ak ' j iW a + W°,i^ k + 
2^%/^ - 2^°^ k \i + (<P akJ + V aj ' k )<P a *° 
-2^ k {^\ aV l a )., - rf k C + f{T ja y k a - f ka ^ a ). (91) 

From (91) and (30), using the relation between y> M „ and in (9), and taking 
into account (34), we find 

int 1__ 1 _ _ _ _ _ _ 

° 3k + t 3k = — - -h J h * - {h?\ ai ¥ k + h k °. m W) 

-(h ai ' j h k a}i + h ai ' k h j a}i ) + (h a3 ' k l + h^iWa + 2h}\ ai h? k + 2h i(r , a h kj }i 
-2h ia , a h H ti + (h ak ' j + h aj ' k )h a<7 '° - 2^ k {h a \Ji\),i - rf k C] 

+ \{T ka h a i + T^h a k ), (92) 

It is easy to verify that total energy- momentum tensor consisting of (17) and 
(92) is conserved. The canonical part of MTW tensor in Hilbert gauge has the 
form 

T 75 = 32^{^' tf V' 7 " l^h' 1 

-rf 5 [\K v ,xh^ x - \h' x h,x]} = T~' 5 - l -^ 5 f. (93a) 
For Newtonian centers this expression coincides with Thirring's one, see eq.(6). 

s 

As noted earlier the nonlocal part of T ^ v is zero for Newtonian centers. 
For these centers from (92) and (9), (12) we have 

e°°+ T 00 = ^(V0) 2 + 2^0, jU = T 00 . (94) 
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For this system the MTW Lagrangian (93) coincide with Thirring's one. The 
same is true for canonical energy-momentum tensors, see (6) and (86-87), but 
spin parts arc different. It follows from (91) and (87) that for Newtonian centers 
MTW spin part contributes only to interaction tensor, while Thirring's spin part 
contributes also to pure field part, see (15). 

We note now that in Hilbert gauge for static case (for Newtonian centers) 
the linearized Einstein equation can be written as 

V 2 C - -V 2 ^, y = -WnGf^, (95) 

see (A12) below. As seen from (93a) for this system T 00 — 0, i.e. there is no 
contribution to hoo from diagrams of Fig. 2. 

Comparing MTW and Thirring tensors in Newtonian approximation we see 
that addition of divergence terms to Lagrangian leads to the change in subdi- 
vision of energy density between purely field part and interaction part. Using 
(61-63) we obtain from (94) the effective gravitational energy density in New- 
tonian approximation given in (56). 

Returning now to contributions from diagrams of Fig. 2 we find that diagram 
(2a) leads to 

ds 2 = (1 + 2<jy)dt -(l-2<f>- <j) 2 )(dx) 2 - tf^-dx'dxK (96) 
Contribution from all 3 diagrams of Fig. 2 is represented in the expression 

i I 

ds 2 = (1 + 24>)dt 2 - (1 - 2<}> - h(j) 2 ){dx) 2 - ^^-dx'dx 3 . (97) 

6 Discussion 

We have assumed in this paper that in 3-graviton vertex each graviton interacts 
with gravitational energy-momentum tensor formed from two other gravitons. 
But in general relativity the 3-graviton vertex is given by cubic in h^ v terms 
in function G(x), see cq. (3), where only quadratic terms are written down. 
Correcting a misprint in [20] one finds that these terms are given by Q^g: 

+h a p n h ali £ — 2h~ fa ^h e l3 ' a + 2h 7a ^h e a '^ + (h ia ^ + hg an )h' a 
—2hj0 ta h' a + (h ja ' a h^ + hga a h^) — h a hfi + 2h 1 e,a.h af3 $ 

^(h^h^^ + hg^h^^)}. (98) 

There are no reasons to expect that Q^g is conserved energy-momentum 
tensor. Moreover it does not contain second derivatives and we know that the 
only conserved energy-momentum tensor with this property is Landau-Lifshitz 
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tensor [16]. Yet (98) does not coincide with LL tensor. This is confirmed also by 
the fact that (98) leads to general relativity result (73), while LL tensor leads 
to (63a). 

We note now somewhat unexpected fact: the half sum of LL and MTW 
tensors reproduces general relativity result (76). On the other hand the half 
sum of Thirring tensor and MTW tensor gives 



In such space-time a relativistic particle moves differently (in G 2 approximation) 
from what is expected according general relativity. 

7 Concluding remarks 

Though general covariance was not assumed, the gauge invariance is of course 
retained [2,10]. For this reason the weak gravitatinal waves in flat space are 
described as in general relativity. All considered tensors give the same energy- 
momentum tensor for the plane gravitational wave. There are no a priori reasons 
to believe that field-theoretical approach will give the same result as general rel- 
ativity. It seems that there is still much to be done to synthesize the geometrical 
and field-theoretical aspects of gravitations. 

I wish to thank V.I.Ritus, and I.V.Tyutin, for useful discussions and D.E.Ivanov 
for collaboration in evaluating the contribution from diagrams (2b), (2c). 

8 Appendix 

Using r/^ — diag(—l, 1, 1, 1) we give here some details of calculating h^ v . We 
utilize the expression 



ds 2 = (l + 2(f>+ 2<t> 2 ) - (1 - 20 + 2<t> 2 ) + 20 




dx l dx 3 



(99) 





where graviton propagator 




r \ 1 f d3 P ■/—> Ou 



(A2) 




P, 



{AA) 



The polarization factor 



satisfies the relations 



1 



P 




p,Vp<7 



tpa ^Vpct — t, 



pLVpO 
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J.IXV p rppo_,llVrp _ \frp _ ,jJ,V rp — T (AK\ 

L r iivpa ± — <- ->- pv 2 — 

The scalar factor D+(a;) has the representation 



and possesses the property 



dr _ 1 

(2n) 2 J_ 00 (x - x') 2 - t 2 + ie 4tt\x - x' | 



drD + (f-5' ,r) = — / - ; - = — - . (At) 



For spherically symmetric body we have to deal with integrals of the kind 

/ d 4 x'D + (x-x')p(r') = -L / = =p(r-Q 

7 47ri y /x>2+x 2 -2x-x> HK ; 



= - / dr'r' 2 p{r') + / dr'r'p{r'), 

r JO Jr 



(AS) 



By the way it is seen from here that the derivative of Newtonian potential over 
r is determined only by the mass inside sphere of radius r. Assuming in (A8) 
that p = ft, we get 

\ jf dr'r' 2 p{r') + jf° dr'r'p(r') = c (± - ^) . (A9) 

Hence the divergent part at small r' appears only in the term, which is absorbed 
by mass renormalization. 

So the source (13) generate the field 



= WnG J d A x'D+{x - x')Too(x') => -<f) 2 . (A10) 



The arrow shows that the divergent part is included in mass renormalization. 
The r.h.s. of (A10) can be obtained also from the solution of wave equation 
derived from (A10) by action of the operator d 2 = V 2 — and taking into 
account that 

-d 2 D+(x-x')=S(x-x') , V 2 ^ = -^. (All) 

We note also that 

V = !6ttG J d 4 x'D + (x - x')f^(x'). (A12) 

Now we shall indicate how to obtain the contribution from diagram (2b). (The 
diagram (2c) contributes as much as (2b) and diagram of Fig.l is irrelevant 
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to finding terms in g^.) First of all we remark that some terms in energy- 
momentum tensor are not symmetric in h forming the tensor. Such terms ought 
to be symmetrized. For example the third term on the r.h.s. of (36) should be 
rewritten as 

-2^, aW => -<^, aW - ip 7 eip af) , af 3- (A13) 

This is important because the left ip will always be considered (for case of 
diagram (2b)) as "contained" in propagator and the right ip as originating from 
Newtonian center. Moreover the second term on the r.h.s. of (A13) may be 
dropped as ip coming from Newtonian center satisfies Hilbert condition (2). 
Rewritten in terms of h (see (7) and (5)) the first term on the r.h.s. of (A13) 
has the form 

-<p af3 \ a p<Pi9 = ~32~£[ /l ° ;/ \aph 7 g - \^h^ a h l6 - ^V-,eh af3 , a ph + ji] 7 gh, a a h] 

(AU) 

Now we consider the contribution from the first term on the r.h.s. of (A14). 
Dropping temporarily the constant factor we have to evaluate the integral 

J d^'lD^x - x')]' a0 h 7 e(x')h^(x'). (A15) 

The last h represents the graviton of Newtonian center. This graviton interacts 
with the source, given by first term on the r.h.s. of (A14). Using (A2), (A4) 
and (69) we come to integrals 

V = - / d 4 x>{D + (x-x%^ 2 (x>), J d 4 x'[D + {x-x% a a <fi 2 {x') = -<fi 2 (x). 

{Al(i) 

The first equation in (All) was used in the last equation in (A16). The first 
integral in (A16) is treated as follows. Integrating by parts we get 

V = j d 4 x'[D+(x - x')}.^')^'). (A17) 

As <fi(x) is independent on i°, v — does not contribute. For the same reason 
fi = also does not contribute. Indeed for \x = we integrate over x'° and get 
the factor 

D + (x-x')\±% = 0, 

see (A6). Thus we replace fi and v by i,j = 1,2,3. Integrating (A17) over x'° 
we obtain 

_ G 2 M 2 f 3 , / 8 1 \ / 9 1\ 1 

h > ~^rj x {dx^W^¥\) 2 {dx^^J 7 {A18 > 

Here we have used (A7) and expression for (fi, see text before eq. (10). Now 
writing 

------- 

r' dx 11 r' ~ 2 dx' 1 r' 2 ' ' 
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and again integrating by parts we find 

G 2 M 2 f 3 , 1_ 8 2 1 

4^rJ \x-x>\dx»dx>i? 2 ~- { ] 

Using relations 

l\ 2Sn 8xV ~{ 2xV 5n\ 

3) --tt + ^ = v2 (-— + -^1 421 ) 

and again twice integrating by parts we get 
G 2 M 2 



7 \ r 4 r- 2 / \x — x'\ 



G 2 M 2 ( -2^- + ^ = - 2r 2 ^^, (A22) 



Finally restoring all factors we find that the contribution from the first term 
in (A14) is 

-4v,0 2 + 8V (A23) 

The divergent integral 



J = j d A x'D + {x-x'){\7 2 cf){x'))cf)(x'), 



appearing in some terms, is cancelled out in final result. 

At last we show how to obtain the finite part of the integral 



(424) 



hi = J d A x'D + {x-x')<P. l3 {x')<P{x') = G 2 M 2 J ^ * 
Utilizing relation 

V 2 I ij {x) = rHlJ r ~~~ C?M\ 
and (cf.(42)) 

2 f ^ X X^ \ X^ 

we find that the finite part of 1^ is 

\t*Mj - S f^ (A25) 

The essential part of this Appendix machinery is checked up by obtaining the 
expression (73) starting from (98). Calculation of 500 can be made by much 
easier method suggested by Schwinger, see [10] and [21]. This method uses more 
fully Hilbert condition and it is helpful for controlling some of our calculation. 
For example it is clear from this method that (A13) do not contribute to h 00 - 
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Figure captions 



Fig. 1 : The second rank tensor formed from matter energy- momentum tensor 
and graviton is a source for another graviton. 

Fig. 2: 3-graviton vertex. Short straight line serves only to distinguish the 
roles of participating gravitons: energy- momentum tensor is formed from two 
gravitons joining the straight line at its ends, this energy-momentum tensor 
serves as a source for graviton emerging from the middle of the straight line. 
Crosses represent external field sources. 
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